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Abstract 

Calculations that demonstrate the influence of three key nuclear reaction rates on the evolution 
of Asymptotic Giant Branch stars have been carried out. We study the case of a star with an initial 
mass of 2Mq and a metallicity of Z = 0.01, somewhat less than the solar metallicity. The dredge- 
up of nuclear processed material from the interior of the star, and the yield predictions for carbon, 
are sensitive to the rate of the 14 N(p, 7) 15 and triple-a reactions. These reactions dominate the 
H- and He-burning shells of stars in this late evolutionary phase. Published uncertainty estimates 
for each of these two rates propagated through stellar evolution calculations cause uncertainties in 
carbon enrichment and yield predictions of about a factor of two. The other important He-burning 
reaction 12 C(a, 7) 16 0, although associated with the largest uncertainty in our study, does not have 
a significant influence on the abundance evolution compared to other modelling uncertainties. This 
finding remains valid when the entire evolution from the main-sequence to the tip of the AGB is 
considered. We discuss the experimental sources of the rate uncertainties addressed here, and give 
some outlook for future work. 
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I. INTRODUCTION 



Modern computer codes for stellar evolution calculations solve routinely a nuclear reac- 
tion network sufficiently large to account for all relevant nuclear transmutations. Reliable 
nuclear reaction rates are a crucial ingredient for accurate modelling of the evolution of 
stars, especially for investigations of the chemical evolution of the stars and of the amount 
of processed material that is returned to the interstellar medium. Such stellar model results 
are used for integrated models of the galactic chemical evolution, and for comparison to 
individual stellar abundance observations. 

In spite of the obvious necessity to check the sensitivity of stellar chemical evolution 
predictions to uncertainties in the underlying reaction rates, little work in this direction has 
been done recently Q|. This is, in particular, true for the evolution of low- and intermediate 
mass stars (0.8 < M*/ M < 8) which host important nuclear production sites and con- 
tribute significantly to galactic chemical evolution. Two main issues make investigations of 
the propagation of nuclear reaction rate uncertainties difficult. One is the computationally 
expensive and numerically difficult nature of the advanced evolutionary phases of low-and in- 
termediate mass stars. The second is that such studies are only useful if consistent estimates 
of the individual rate uncertainties are available. 

We have started a program to address this problem. In Paper I we investigated the 
impact of CNO cycle (p, 7) reaction rate uncertainties on the predicted stellar oxygen isotopic 
ratios. These predictions are important for the astrophysical interpretation of pre-solar 
meteoritic corundum grains P,Q|. In Paper I we used a Monte Carlo approach integrated 
into nuclear network post-processing calculations. 

In this second study we focus on the evolution of Asymptotic Giant Branch (AGB) stars 
and their sensitivity to the rates of three key nuclear reactions: 14 N(p,7) 15 0, triple-a and 
12 C(a, 7) 16 0. We extend our preliminary results of this project ^], and provide a more 
in-depth presentation. The following sections describe: Sect. |H] astrophysical background; 
Sect . lIHI nuclear physics input, its uncertainties and possible revisions; Sect.HYlthe physical 
model and methods of the astrophysics simulation, as well as the main elements of AGB 
evolution that are important here; Sect.|Sour results as well as additional calculations with 
a second, independent code that verify the findings; Sect.|Vl] new results for the triple-a 
rate; and Sect. lVIll results and discussion. 
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II. ASTROPHYSICAL CONTEXT 



After the initial H- and He-core burning; phases, low- and intermediate mass stars evolve 
into dooMe-shell-burning giant 8tarS fl^ These AGB 8tarS have a u„ lqU e m echanical 
structure, as explained in more detaniaHe^B.Theelect.on-degen.ate core with mass 
M c « 0.6 M consists of carbon and oxygen (the ashes from previous evolutionary phases) 
and is about the size of the earth, while the envelope has on average roughly the density of 
water and extends to several hundred times the radius of the sun. This envelope is unstable 
against convection, and is thus well mixed. 

Recurrent He-shell flashes, with periods of 5 to 10 xl0 4 yr are a characteristic of these 
configurations. A combination of partial degeneracy, the small geometric scale of the He- 
burning shell, and the strong temperature dependence of the triple-a reaction rate leads to a 
thermonuclear runaway that locally generates power of roughly 10 34 watts (10 8 L ). Neither 
heat conduction nor photon radiation are sufficient to carry away the energy, and the layer 
between the He- and the H-shell becomes convectively unstable. 

Part of the energy released in the He-flash will do expansion work, cooling the layers above 
the He-shell. The stellar opacities in this region at the base of the convective envelope then 
increase, which forces the envelope deeply into the core (in terms of the Lagrangian mass 
coordinate, the enclosed mass). This penetration of the envelope convection zone into the 
processed core material below the the H-burning shell is the third dredge-up. 1 The details 
of these events are well documented in the astrophysical literature 0, liol . 
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The dredge-up event is responsible for the transfer of nuclear processed material from 
the high-temperature stellar interior to the low-temperature stellar surface where it can 
be observed spectroscopically. This material is also blown into the interstellar medium by 
stellar winds. For these reasons, the strength of the dredge-up is of great importance to 
the observed chemical enrichment of low-mass giants and their role in galactic chemical 
evolution. 

The amount of dredge-up obtained in stellar models of the AGB stars continues to be a 
matter of debate. It is well established that the dredged-up amount depends on the core 
mass, the stellar metallicity and opacity, the model of convection, and the treatment of 



1 It follows two dredge-up periods during previous evolutionary phases, which are however less important 
for the overall chemical enrichment of low- and intermediate mass stars. 
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convective boundaries, as well as on the numerical implementation of several details in the 
codes Earlier studies have already indicated that stronger He-shell flashes are 

followed by deeper dredge-up )l8|. and that a decreased energy generation in the H-shell 
leads to stronger He-shell flashes [19| . However, there are no investigations of the sensitivity 
of dredge-up and the envelope abundance evolution to nuclear reaction rate uncertainties. 

We focus on the three reactions that dominate H-burning and He-burning respectively. It 
is well known that the 14 N(p, 7) 15 reaction has the smallest rate in the CN cycle; it controls 
the circulation rate of the CN catalytic material in the H-burning shell. The 12 C(a,7) 16 
reaction is weakly active during the interpulse phase between the He-shell flashes. The most 
important reaction for the He-shell flash is the triple-a reaction. 

In double-shell burning around degenerate cores several competing time-scales are in- 
volved. The period of the He-shell flashes depends mainly on the rate of Helium accretion 
from the H-burning shell. The strength of the He-shell flash depends on the geometrical size 
and on the partial degeneracy of the He-shell. A more degenerate and thinner shell results 
from a longer flash period, which in turn can be caused by the smaller energy generation 
rate due to a smaller CNO cycle rate. In such a case it takes longer to accrete the required 
amount of He from the H-burning shell to ignite the flash. Thus, one can qualitatively 
understand that the H-shell burning rate influences the He-shell flash strength, and thereby 
the subsequent dredge-up. 

The possible influence of the triple-a reaction on the He-shell flash strength is perhaps 
more obvious. A larger rate is likely to cause a larger peak-flash He-burning luminosity, 
and subsequently a deeper dredge-up. Our qualitative expectation is then that a reduced 
14 N(p,7) 15 rate and an increased triple-a rate will each increase the amount of carbon 
produced in the process. The effect of the 12 C(a,7) 16 rate is less obvious. To obtain 
quantitative estimates on these processes we have conducted a detailed numerical study. 



III. NUCLEAR PHYSICS INPUT 

n 

The NACRE collaboration |20| has recommended reaction rates for the reactions we 
consider here. In this section we examine whether, six years after their publication, the 
NACRE estimates still describe the available data with sufficient accuracy for our purposes. 
Table H] lists the temperatures of interest: T§ = 0.5 (T = 5 x 10 7 K) to T§ = 0.8 for the 
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14 N(p,7) 15 reaction and T 8 = 1-3 for the triple-a and the 12 C(a,7) 16 reactions. 
A. 14 N(p, 7 ) 15 

The 14 N(p, 7) 15 reaction has a complex structure with transitions, both resonant and 
non-resonant, to several final states contributing to the rate. In addition, contributions from 
the tails of subthreshold states must be considered. Schroeder et al. |2l( obtained data over 
a wide energy range from 0.2 to 3.6 MeV; their R-matrix fit to this data yielded a large S- 
factor for the transitions to the ground state. Their total S-factor, S(0 ) = 3.2 ± 0.54 keVb, 
was the principal basis for the NACRE reaction rate 3.2±0.8keVb 20j. However, the fit to 
the ground state cross section required an unusually large value for the gamma width of a 
subthreshold 3/2 + state at E x = 6.793 MeV in 15 0, about 7 times that of the isospin-analog 
ransition in 15 N. Such large differences are seldom, if ever, seen, at least for light nuclei 



22, 23. . Motivated by this fact, direct measurements of the 6.793 state's lifetime were made 
2J, |25j and yielded much smaller gamma widths. A reanalysis [26] of the Schroeder data 
resulted in a much smaller ground state transition and 5"(0) = 1.77 ± 0.2 keVb. Later, 
determinations of asymptotic normalization coefficients (ANC) for the relevant transitions 
using nuclear transfer reactions, combined with some of the Schroeder results led to 
S(0) = 1.7 ±0.41 keVb. 

Recently, additional 14 N(p, 7) 15 data and corrected data from Schroeder et al. [2]] were 

28| . An independent measurement at TUNL 



analyzed to yield S(0) = 1.70 ± 0.22 keVb 

(2^] yielded S'(O) = 1.68 ±0.18 keVb. Measurements of the analyzing power at 270 keV 
indicate that Ml contributions to the cross section should be considered; to our knowledge 
these contributions have not been included in detailed fits to the data. There are also 
preliminary data down to 70 keV 31] for the main transition, that to the 6.79 MeV state. 

All the recent investigations show that the ground state transition is small and that the 
resulting total S factor is smaller by about a factor of two than the NACRE result. In 
the near future, the reliability of the rate is likely to improve as more complete analyses 
including all the recent data are carried out. However, for now we have chosen to use an 
unweighted average of the four recent results and a conservative error reflecting the relatively 
long extrapolations to the astrophysical range, the neglect of the Ml amplitudes, and the 
imperfect fits to the data. We obtain S(0) = 1.70 ± 0.25 keV b. For easy employment in 
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the stellar evolution code, this value can be approximated as a fraction / = 0.64 ± 0.1 of 

n 

NACRE's analytical fit [20j to the reaction rate in the relevant temperature range. 



B. Triple-a 

The first step of the triple-a process is the fusion of two a particles to form a equilibrium 
concentration of 8 Be. The subsequent capture of an alpha particle produces an equilibrium 
concentration of 12 C in its 7.65 MeV + state. Occasionally 12 C is formed by a leak, via a 
gamma cascade or pair emission, to the ground state of 12 C. For the temperatures involved 
in the present calculations both of these steps are resonant and the reaction rate is given by: 



r^a oc T rad exp(-Q/kT). (1) 

The value of Q for the 7.65 MeV state is known to within ±0.2 keV and contributes 
an uncertainty in the rate of only ±1.2% for T§ = 2. Essentially all the uncertainty in the 
rate is due to the uncertainty in T ra< i, the radiative width of the 7.65 MeV state; T ra d is 
known with a precision of ±12%. It is essentially these established values of r ra d and Q that 
are incorporated into the NACRE rates for the temperatures considered here, and hence, 
for our purpose the NACRE rates are adequate. We shall see, however, that they are not 
sufficiently accurate. 

Following the completion of the calculations described herein, H. O. U. Fynbo, et al. j^, 
determined the level structure of 12 C and concluded from their results that the 7.65 
MeV state alone adequately describes the reaction rate for temperatures T = 0.1 — 100 x 
10 8 K. They also conclude that at T = 2 x 10 8 K, the midpoint of the temperature range 
covered here, the triple-alpha rate is smaller than the NACRE results by about 10%. This 
is comparable to the quoted error we used, and should be kept in mind when examining the 
details of the present results. In fact, our calculations can be used to determine the effect 
of the 10% change on C-yields. 



C. 12 C(a, 7 ) 16 

This reaction has been the subject of many experiments and analyses over a period of 
forty years, but it is still not accurately known. It is not possible here to review this subject 
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in detail, it is simply too complex. One can find a comparison of the various rates in Fig. 1 
of Eid et al. [35], and the references cited there can be consulted for more detail. A extensive 
discussion of this rate can be found in Buchmann 36( and the results of extensive recent 



measurements in J. W. Hammer, et al. |37|. The NACRE rate is probably somewhat too 
large in the region of present interest, but the quoted uncertainties are sufficiently large 
to represent the probable range of acceptable values. This is a minor issue for the present 
calculations since, as we shall see, carbon production in low and intermediate mass AGB 
stars is very weakly dependent on the 12 C(a,7) 16 rate. 

D. Comments 

We conclude that the NACRE recommendations describe current triple-a and 
12 C(a,7) 16 data sufficiently well, at least for the temperature range relevant here. Re- 
cent data, however, make it clear that the NACRE estimate for the 14 N(p, 7) 15 rate is high 
by roughly a factor of two. 

IV. METHODS AND PHYSICS INPUT 

The one-dimensional stellar evolution codes we employed solve the well-established full 
set of non-linear partial differential equations to account for hydrostatic equilibrium, mass 



14. 



16 



3811 . Most of the calculations have 



continuity, and energy transport and generation 

been done with the code EVOL jsjj. It includes updated input physics. 2 A small amount 
of envelope overshooting, but no overshooting at other AGB convection zones is introduced 
following Herwig et al. 42]. 3 

In order to generate the initial model for our comparative study we compute the evolution 
of a star with a mass of 2M and a metallicity of Z = 0.01 from the pre-main sequence 
through the H- and He-core burning phase. After the He-core burning phase the star grad- 
ually climbs up the Asymptotic Giant Branch in the logT e ff - logL (Hertzsprung- Russell) 



2 The opacities, for example, are from Iglesias and Rogers |4p| supplemented with low temperature opacities 
by Alexander and Ferguson 01- 



3 The efficiency for convective envelope (CE) overshooting is f cc = 0.016, see Herwig for a description 
of the overshooting scheme used. Mass loss is included by adopting the formalism of Blocker with a 



scaling factor ijb — 0.1. For more details, and definitions see Herwig [39j. 
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TABLE I: Relevant temperature range, NACRE nuclear reaction rates and their uncertainties 
20], and adopted factors to fitting formula rates for our calculations. 



reaction T 8 < av >i ow < av > a < av > h i g h exp 6 f fit c f up d f iow e 

14 N(p, 7 ) 0.5 2.67 3.68 4.69 -10 0.9701 1.3137 0.7479 

0.8 0.76 1.04 1.32 -7 0.9488 1.3377 0.7702 

adopted for our calculations: 1.33 0.75 

3a 1.0 2.05 2.38 2.70 -24 1.0424 1.0883 0.8263 

3.0 3.95 4.57 5.18 -13 1.0068 1.1258 0.8585 

adopted for our calculations: 1.13 0.82 

12 C(a,7) 1.0 1.06 1.81 2.55 -20 0.9762 1.4431 0.5999 

3.0 2.88 4.75 6.62 -12 0.9905 1.4070 0.6121 

adopted for our calculations: 1.44 0.60 



"Recommended reaction rate. 

kpower of 10 multiplying reaction rates in columns 3,4, and 5 
c Ratio between tabulated value and fit formula. 

d < av > hig h /(< av > ff it ). 
e < av >iow /(< av > ffu) 

diagram. We choose as a starting model for all subsequent calculations a model at the very 
end of the He-core burning phase, and well before the onset of the first He-shell flashes. 

The evolution of the H- and He-burning shell and of the stellar surface in the Lagrangian 
mass coordinate for the benchmark sequence ET2 are shown in Fig.^ From the end of the 
He-core burning phase the star spends about 20 million years on the so-called Early-AGB 
phase. During this phase the H-burning shell is largely inactive and most nuclear energy is 
produced in the He-burning shell. 

He-shell flashes occur only during a rather short period of the post-He-core burning phase. 
The underlying reason for their occurrence is the different burning rate of the two shells which 
eventually prohibits quiescent double-shell burning. A close-up of the actual He-shell flash 
phase of the AGB is shown in Fig.|21 Seventeen thermal pulses can be identified by the 
vertical lines that connect the H- and the He-shell at almost equidistant intervals. These 
vertical lines are the brief He-shell flash convection zones which last for only 200 to 300 yr. 
During the flash the convectively unstable layers are confined to the region below the H-shell 
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and above the He-shell. The inset shows a small spike at the bottom which is the rapidly 
growing upper boundary of the He-shell flash convection zone. It stops just short of the mass 
coordinate of the H-free core where the H-shell is located, and quickly retreats afterward. 
The H- and He-burning shells remain well separated even during the He-shell flash episodes, 
and no H from the envelope can enter the He-burning shell. The inset also shows how the 
bottom of the convective envelope later descends into mass layers previously occupied by 
the He-shell flash convection zone and with consequent "dredging" of processed material 
into the envelope. It is this tiny detail in the convective evolution of the stellar interior that 
is responsible for the enrichment of the envelope and eventually, through mass loss, of the 
interstellar medium. The dredge-up events after the thermal pulses cause a gradual increase 
of carbon and to a much lesser extent oxygen (Fig.|3J). The surface abundance for 16 is 
nearly constant with time for near-solar metallicity. This results in an increase of the C/O 
ratio, and eventually to the formation of C-stars. More detailed figures of the evolution of 
He-shell flashes can be found, for example, in Fig. 1 and 10 in Herwig 

For any comparative study of the propagation of nuclear reaction rate uncertainties in a 
stellar evolution code a somewhat consistent set of quantitative estimates on the uncertain- 
ties is required. In the absence of such estimates, one is left with the rather crude approach 
of applying common factors to all rates in question. However, such an approach misses out 
on critical aspects of the error propagation in a real stellar evolution environment. As we 
will show, for example, the reaction in our sample with the largest relative error has the 
smallest impact on the observable prediction. Similar conclusions were drawn in Paper I Q] . 

In this study we rely initially on the NACRE compilation 20] which contains recom- 
mended values and estimates for lower and upper bounds as a function of temperature. We 
then consider the impact of the revised recommendation for the 14 N(p, 7) 15 reaction. Note 
that the recommended factors only apply in the given temperature range. In Tabled and 
Fig.0] we show the relevant information from the NACRE compilation for the two temper- 
ature ranges appropriate for He- and H-burning respectively. In addition to the tabulated 
reaction rates, fitting formula for the recommended values are provided. In the calculations 
we use these formulae instead of tables to evaluate the reaction rates at the required tem- 
peratures. We have checked the accuracy of the formulae, and found that for the T range 
of interest here (as indicated in Table|lJ column 2) the fitting error is rather small as shown 
in Column 7 of Tabled which shows / fit =< av > ta bie / < crv > fit formula ■ 
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TABLE II: Results for NACRE recommended rates and uncertainties. 



ID 


reaction 


factor 


iV T P a 


X(m c = 0.56 M 


\ b \ c 

J A max 


EM dup /lO- 2 M d 


pi2 C e 


ET2 


all 


1.00 


8 


0.16 


0.29 


1.2 


2.19 


ET5 


14 N(p, 7 ) 


1.33 


9 


0.15 


0.31 


1.2 


1.95 


ET8 


14 N(p, 7 ) 


0.75 


11 


0.29 


0.39 


2.3 


4.27 


ET6 


3a 


1.13 


10 


0.31 


0.41 


2.5 


5.42 


ET9 


3a 


0.82 


7 


0.12 


0.29 


1.1 


1.79 


ET7 


12 C(a, 7 ) 


1.44 


9 


0.24 


0.34 


1.6 


2.71 


ET10 


12 C(a, 7 ) 


0.60 


8 


0.21 


0.33 


1.5 


3.12 


"Number of Thermal Pulses that cause dredge-up. 

^Dredge-up efficiency A at mass coordinate 0.56 M , for details see text. 
c Maximum A reached by any flash in the entire sequence. 
d Entire mass dredged-up by all dredge-up events. 

el2 C yield from thermal pulses and the accompanying dredge-up in units of M Q . For details see text. 






TABLE III: Results for revised 


14 N(p, 7 ) 15 


range. 




ID 


14 N(p, 7 ) 


3a 


N TP a 


A(m c = 0.56 M 


) ^max 


£M dup /lO- 2 M 




ET12 


0.75 


1.13 


10 


0.33 


0.43 


2.5 


5.65 


ET13 


0.64 


1.00 


9 


0.29 


0.41 


2.2 


4.62 


ET14 


0.64 


1.13 


11 


0.36 


0.44 


2.9 


6.02 


ET15 


0.54 


1.13 


11 


0.41 


0.46 


3.2 


6.98 


ET17 


0.54 


0.82 


9 


0.25 


0.39 


2.3 


5.29 


ET18 


0.75 


0.82 


8 


0.16 


0.41 


1.8 


3.98 



a See TableTm and text for explanations and details. 

V. CALCULATIONS AND RESULTS 

A. Chemical enrichment and dredge-up as a function of nuclear physics input 

Starting from our initial model at the end of core He-burning, we calculate seven full 
evolutionary sequences which end when all envelope mass is lost and the remaining stellar 
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core is about to become the central star of the planetary nebulae stage. The sequences 
differ only in the adopted rates for the three reactions investigated here fTable lTT|) . One 
benchmark sequence (ET2) has been computed with the recommended NACRE rates for all 
three reactions. In addition six sequences have been calculated in which for each of these 
reactions the rate from the fitting formula is multiplied by the factors given in Tabled which 
include the small differences between tabulated and fitting formula values. Thus, for each 
reaction two sequences are computed, one adopting the upper and one adopting the lower 
bound of the uncertainty range. 

In a second set of sequences we have investigated the influence of the revised (smaller) 
14 N(p,7) 15 rate. We have considered also simultaneous changes of both this rate and the 
triple-a rate fTable lTTTJ) . An overview of the reaction rate choices for the most interesting 
cases in which the 14 N(p, 7) 15 and the triple-a rate have been changed is shown in Fig.0] 

We summarize these results in TablelTTl and Table HTTl The most salient features of these 
results are described by the efficiency of dredge- up, A, and the yield of carbon pi2 C . The 
efficiency is given by A = AMdup/AMn with AMh the core mass growth between two He- 
shell flashes due to H-shell burning, and AMd up the dredged-up mass following the He-shell 
flash. For each flash A indicates the efficiency of dredge-up. A = 1 means that the same 
amount of mass by which the core grew between two flashes is dredged-up after a flash. 
The 12 C yield from thermal pulses and the accompanying dredge-up is given for i = 12 C by: 
Pi = /^(Xj(m) — Xi n i) dm , where M\ and Mf are the initial and the final stellar mass at the 
beginning and the end of the AGB phase, X is the mass fraction at the surface as the star 
evolves and X in \ is the initial mass fraction. We also tabulate the total dredged-up mass. 

The chemical evolution of AGB giants depends sensitively on the rates of the 14 N(p, 7) 15 
and the triple-a reaction. Calculations with a smaller 14 N(p,7) 15 rate (case ET8) show a 
larger dredge-up efficiency, a larger dredged-up mass, and a larger amount of carbon mixed 
from the processed layers to the envelope and a larger carbon yield. 4 All these quantities 
are about a factor of two higher than for the benchmark case. 

For the triple-a reaction we observe the opposite behavior. The case ET6 with a larger 
rate has, on average, He-shell flash peak luminosities (not shown in the table) that are 20- 

4 The carbon enrichment is used here as a proxy for the envelope enrichment with nuclear processed material 
which would include the s-process elements, for example. 



12 



50% higher than the benchmark case and accordingly dredge-up is deeper. The efficiency, 
total dredged- up mass, and the 12 C yield are about a factor of two higher for sequence ET6 
compared to sequence ET2. It is also clear that these effects are non-linear in the rates. The 
changes for increases in the 14 N(p,7) 15 rate and decreases in the triple-a rate are small. 

The uncertainty in the 12 C(a, 7) 16 as given in the NACRE compilation rate has a much 
smaller influence on the observables studied here. We find that in both cases with upper 
and lower range values for this rate (cases ET7 and ET10) the 12 C yield is somewhat larger 
than in the benchmark case. Sequence ET10 has a slightly larger 12 C yield than ET7 
although ET7 has a slightly larger dredge-up mass and efficiency than ET10 (see discussion 
in Sect. EE) . 

In Table lTTTI we show the results for the revised 14 N(p, 7) 15 rate and for cases in which two 
rates are changed simultaneously. The 12 C yield increases further for the lower 14 N(p,7) 15 
rate and the combination with an increased triple-a rate leads to still larger values. A 
compact representation of the set of simulations is given in Fig. El The main result is that 
due to the revision of the 14 N(p,7) rate the predicted 12 C yields of low-mass stars is about 
twice as large as with the old rate. In addition, the relative error of the combined effect of 
the 1 N(p, 7) and triple-a rate has decreased. 

The evolution of the 12 C abundance is useful to further illustrate the differences. In Fig-EI 
the envelope 12 C abundance increases in discrete steps for all cases. These steps correspond 
to the discrete dredge-up events after sufficiently strong He-shell flashes. The astrophysical 
yield is obtained by integrating the surface abundance over the mass lost. Most notable 
is the fact that 12 C abundances for the shown cases span a range significantly exceeding a 
factor 2. Until the reaction rates are better known this is an unavoidable uncertainty in the 
yield predictions. 

These results confirm our original qualitative expectations. A smaller 14 N(p, 7) 15 rate 
leads to a smaller helium production rate and later ignition of the He-shell flash. This flash 
is then more violent and the subsequent dredge-up is more efficient compared to a case with 
a larger 14 N(p, 7) 15 rate. More efficient dredge-up leads to a larger envelope enrichment, 
and thus the envelope 12 C abundance and stellar yield is larger for the run with the smaller 
14 N(p, 7) 15 rate. For the triple-a reaction a larger rate leads to stronger He-shell flashes. In 
fact the run with the large rate shows He-burning peak luminosities which are about 20-50% 
larger than for the run with the lower triple-a rate. Accordingly the run with the larger rate 
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shows greater efficiency, deeper dredge-up, and larger 12 C abundances at the surface and in 
the yields. 

B. Verification of stellar modelling results 



The main set of models has been computed using the EVOL stellar evolution code [16j . In 
order to check these results we have repeated a subset of the numerical experiments with the 
independent MSSSP code jllj]. Calculations were carried out for M = 2.1 M Q and Z = 0.008 
with the standard 14 N(p, 7) rate used in that code (very similar to the NACRE rate) starting 
from the main-sequence and a comparison calculation starting, as with the EVOL set of 
sequences, after the end of core He-burning with 0.6 times the standard 14 N(p, 7) rate. This 
calculations shows deeper dredge-up and a larger C yield compared to the sequence with 
the higher 14 N(p,7) rate. The effect seen in the MSSSP calculations is qualitatively and 
quantitatively consistent with the EVOL results. 

In several test calculations we discovered that starting the comparison runs with different 
reaction rates at or even after the onset of the He-shell flashes did not show consistent trends. 
The luminosity of the H- and He-shell depend on the core properties (im particular mass 
and radius) that are the result of the previous core burning phases of H and He. For the 
main set of comparative calculations we therfore choose an initial model about 20 million 
years before the first thermal pulse, and immediately after the end of the He-core burning. 

Comparison calculations have been done with the MSSSP code with different 12 C(a, 7) 16 
rates. As with the EVOL calculations reported in Sect . IV Al we were not able to identify 
clear correlations between this reaction rate and the dredge-up and yield properties of the 
models. The differences we found were below the 10 — 15% level and more sensitive to 
numerical parameters (such as the spatial and temporal resolution) than the 12 C(a,7) 16 
rate. (This is actually to be expected, because this reaction rate is not as important as the 
triple-a rate during the AGB thermal pulse evolution.) The quantities like total dredge-up 
mass and carbon yield are the result of a discrete process repeated n » 8 to 11 times. Thus, 
these numbers are subject to a statistical fluctuation ~ 1/n. 
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C. The pre-AGB evolution: H- and He-core burning 

The main emphasis of this study is the impact of nuclear reaction rates on the AGB 
evolution, in particular the chemical enrichment through the third dredge-up. However, we 
did check how the dredge-up of models with different reaction rates change if run all the way 
from the main-sequence, including the effect on H- and He-core burning. One test has been 
done with the MSSSP code, starting a sequence with the reduced 14 N(p,7) rate from the 
main-sequence. In that sequence the dredge-up during the AGB thermal pulse phase is still 
significantly larger than with the NACRE recommended rate, but the increase is somewhat 
smaller then for the comparison runs started after the end of He-core burning. 

A second test has been made with the EVOL code running two sequences all the way 
from the zero- age main-sequence to the end of the thermal pulse AGB. In addition to the 
benchmark case ET2 we reran the combination of case ET14 (0.64 x 14 N(p, 7), 1.13 x 3a), 
and refer to this run as ET14a. We found that in this case the dredge-up is about 20% 
larger than in the ET14 model calculated from the starting model after the end of He-core 
burning. This is another example of the highly non-linear behaviour of the third dredge-up. 

We can analyse the differences caused by the change of nuclear reaction rate on the 
H- and He-core burning phases, and find them to be very small. The central temperature 
during the H-core burning of the 2 M Q models studies here are 2.1 x 10 7 K initially, increasing 
sharply to 3.45 x 10 7 K at the end of H-core burning. Case ET14a shows central temperatures 
throughout core H-burning that are 1% larger than the benchmark case. This small increase 
of temperature is sufficient to increase energy generation required for hydrostatic equilibrium 
because of the steep temperature dependence of the 14 N(p, 7) rate. Run ET14a consumes H 
in the center slightly faster, and accordingly the H-core burning phase is about 1% shorter 
than in the benchmark case. The mass of the convective core is practically the same in both 
cases. None of these nuclear reaction rate differences during the main-seqeunce evolution 
would make an observable difference. 

During the He-core burning the slightly larger triple-a rate leads to a slightly larger C/O 
ratio on the core (2%). However, during the EVOL calculations some breathing pulses of the 
convective core occur. These breathing pulses are well known during the He-core burning 
phase and related to the unstable growth of the convective core into a layer that is stabilized 
by an composition gradient. The treatment of convective boundaries implemented in the 
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EVOL code makes the occurrence of breathing pulses somewhat dependent on the numerics. 
It is not clear from the calculations to which extent the magnitude of the central C/O ratio 
depends on this simulation error. However, since the effect is small we decided not to follow 
this question any further. At the end of He-core burning the age difference between the two 
runs is 0.3% of the total age. The core mass and size are practically identical. 

The second dredge-up decreases the core mass slightly for stars of this mass. This effect 
is weaker in run ET14a, so that the core mass after the second dredge-up is slightly larger 
than in the benchmark run. It is during the early AGB evolution from the end of He-core 
burning to the first thermal pulse that sequence ET14a has a slower growth rate of the H- 
free core, corresponding to a 3% smaller H-burning luminosity. This leads to a 2% smaller 
core mass at the first thermal pulse for case ET14a compared to the benchmark run. A 
smaller core mass should result in less efficient third dredge-up, if the reaction rates are the 
same. However, ET14a has a combination of rates that increase third dredge-up, as shown 
in the main set of comparative calculations. In addition, because the stellar luminosity of 
run ET14a is somewhat smaller, the mass loss according to the adopted L-dependent mass 
loss formula is smaller, and the thermal pulse AGB phase is longer by 15%. In particular 
the last thermal pulses have larger He-flash peak luminosity, that as discussed above leads 
to more efficient dredge-up. 

The third test concerns the influence of the 12 C(o;,7) 16 rate on the pre-AGB evolution. 
We calculated an additional sequence (ET7a) corresponding to sequence ET7 (NACRE rate 
times 1.44) from the main-sequence to the thermal pulse AGB phase. As can be expected, 
there are no difference (e.g. duration, core size) between run ET7a and ET2a during the H- 
core burning phase. During the core He-burning phase the central C/O ratio is systematically 
smaller in run ET7a compared to the benchmark case. For example after 1/3 of the He-core 
burning phase C/O = 3.6 for the ET7a run and C/O = 5.2 for the benchmark case. Not 
surprisingly the ratio of those two values is 1.44. The ET2a He-core burning duration is 
3% shorter than ET7a, the central temperature and density are almost the same. As in the 
discussion of run ET14a above we note the occurence of breathing pulses during the He-core 
burning phase, in particular towards the end of this phase. In runs ET2a and ET7a these 
breathing pulses are very similar, with no observable difference. The core mass at the first 
thermal pulse is 0.494 M for ET7a, and practically the same (0.493 M ) for ET2a. We 
continued run ET7a into the thermal pulse regime until the onset of the third dredge-up. 
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As expected from the practically idenitical core masses the ET7a and ET2a sequence show 
very similar third dredge-up behaviour. 

In conclusion we find that as expected the sensitivity of stellar evolution properties to 
the 12 C(a, 7) 16 rate during the H- and He-core burning phases of low-mass stars are small. 
There is a small dependence of the thermal pulse AGB results on the progenitor evolution, 
which however does not change the trends established with our main set of comparative 
calculations starting with the same initial model after the end of He-core burning. 



VI. IMPROVING THE TRIPLE-a RATE 

The recent results of H. O. U. Fynbo, et al. (^j], j3| have shown that for T = 0.1 — 100 x 
10 8 K, the triple alpha reaction rate is dependent essentially on the properties of the + 
state at an excitation energy 7.65 MeV in 12 C, the Hoyle state. This reduces the problem 
of determining the triple alpha rate to determining the properties of that state. 

Because the 7.65 MeV state is a + state its direct excitation is difficult. Moreover, the 
ratio of the radiative width to the total width is small, 4.13 x 10~ 4 . As a result one must 
determine T ra d from the relationship 

r + r r 

r — r _i_ r — 7 _ r (o\ 

1 rad — 1 7 ~r J- n — p r ^ ' 

I 1 

Here T 7 , T n and T are the gamma width, pair width and total width of the 7.65 MeV state. 
Each of the three factors on the right is determined in a separate experiment. At present 
they are known with an accuracy, left to right, of ±2.7%, ±9.2%, and ±6.4%. In all cases, 
there are several consistent measurements, so these results can be regarded as robust. 

There are two new developments that may significantly improve our knowledge of T ra d, 
and hence the triple-a reaction rate, by improving the accuracy of the poorest known quan- 
tities: the pair width T v and the pair branch T^/T . The pair width is determined from the 
transition charge density for inelastic electron scattering to the 7.6 MeV state. There is a 
new, as yet unpublished result based on a compendium of extant measurements over 
a large momentum transfer range, that has a quoted accuracy of ±2.7%. It is difficult to 
imagine that a more accurate value of T n can be obtained. On the other hand, this value is 
not quite consistent with the earlier values of T n . 

The pair branch T n /T is the least well known quantity, primarily because it is so small, 
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about 6 x 1CU 6 . A new experiment Q], a Western Michigan University (WMU), Michigan 
State University (MSU) collaboration, is underway using the Tandem accelerator at WMU. 
The proposed detector is an improved version of that used by Robertson et al. 45]. 

In this experiment the 7.6 MeV state in 12 C is excited by inelastic proton scattering, 
taking advantage of a strong resonance at an excitation energy of 10.6 MeV and a scattering 
angle of 135 degrees in the lab. In order to reduce gamma ray backgrounds, a coincidence 
is required between a thin plastic- scintillator cylinder surrounding the target and a large 
plastic scintillator surrounding both the target and the cylinder. This arrangement should 
strongly discriminate against 7 ray-backgrounds-gamma rays have only small probability of 
interacting in the thin cylinder. The pair branch is then given simply by the ratio of the 
number of positron-electron pairs detected by plastic scintillator coincidences to the number 
of counts in the 7.65 MeV peak in the proton spectrum. An examination of the systematic 
uncertainties in the similar Robertson experiment leads us to estimate that an accuracy of 
5% is achievable. 

These two results promise to reduce the uncertainty in the triple-a rate to about 6%; as 
we have seen that will greatly improve the reliability of predictions of carbon production in 
AGB stars. 



VII. DISCUSSION 

We have presented a systematic investigation of the propagation of the rate uncertainties 
of key nuclear reaction rates into chemical enrichment predictions of low- and intermediate 
mass stars that have reached the thermal pulse AGB phase. We found that the dredge-up 
in low-mass stars depends rather sensitively on the adopted reaction rates. The overall 
dredge-up of material and, specifically, the yield of 12 C, has uncertainties of greater than a 
factor of two owing to the reaction rate uncertainties. The C/O ratio at the stellar surface 
has a similar uncertainty. 

Such uncertainties are a problem for many problems of current astrophysical interest. 
The construction of integrated models of galactic chemical evolution, for example, includes 
contributions from stars of all initial masses , and AGB stars are an important contributor 
for some nuclear species. The enrichment of the surface abundance with carbon also effects 
the appearance of AGB stars in extra-galactic stellar population studies. As the surface 
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abundance changes from O-dominated (C/O < 1), to C-rich (C/O > 1) the molecular 



chemistry in the giant's atmosphere changes considerably [47[, affecting the star's surface 
temperature and thereby its astronomical colors. In older extra-galactic populations AGB 
stars are often the brightest stars, and can probe the population's properties, for example, 
its age. Finally, many extremely metal-poor (EMP) stars, which may provide information on 
chemical evolution in the early Universe, turn out to have binary White Dwarf companions 



|48j . The unusual abundance patterns of these EMP stars [49] should correlate with the 
chemical yields of the White Dwarf progenitors - the AGB stars at this low metallicity. 

Of course nuclear reaction rates are not the only uncertainties in AGB models. The 
treatment of convection, and mixing in general, effects the efficiency of the third dredge-up 
as well 0, . Two separate issues have to be considered. Convection in ID stellar evolution 
models is usually approximated by some variant of the local mixing-length theory j^|.In this 
ballistic theory the mean free path of rising and descending blobs has to be specified. 5 Based 
on the sparse information on this topic in the literature Q 0, we roughly estimate 
that the mixing-length uncertainty translates into yield uncertainties ranging from 30% to 
a factor of a few, depending on initial stellar mass. 

Another source of uncertainty of dredge-up predictions is the treatment of convective 
overshooting. There is now enough numerical and experimental proof to claim that convec- 
tive overshooting takes place in stellar environments, and that the efficiency of that process 
depends on the evolutionary phase (3). It appears that dredge-up predictions are uncertain 
by a factor of two because of the poorly known overshooting efficiency. 

These two issues related to the modelling of stellar convection have been viewed as the 
major source of dredge-up and yield prediction uncertainties. Our study shows that nuclear 
reaction rate uncertainties of two key reactions induce modelling uncertainties of similar 
magnitude. The need to reduce these uncertainties is a powerful argument for better de- 
terminations of the reaction rates of the 14 N(p, 7) 15 and triple-a reactions. Progress in 
experimental nuclear physics will have an immediate impact on astrophysical models that 
rely on stellar yields. 



5 Usually one uses the well known stellar parameters of the sun to calibrate this free parameter, and keeps 
this value constant as the evolution progresses. However, multi-D hydro-simulations have shown that this 
assumption is not correct |52( , and for evolved giants the mixing-length parameter may be larger by some 
substantial fraction. 
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FIG. 1: Evolution of Lagrange coordinates of stellar surface, H- and He-burning shell with in- 
creasing time. As the star evolves mass loss increases with time and the total mass decreases. The 
lines start at the end of He-core burning. Time has been set to zero at the maximum He-burning 
luminosity of the first He-shell flash. 
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FIG. 2: Similar to Fig.^ but only for the AGB phase with He-shell flashes. The inset shows the 
detail of an individual dredge-up event. The inset mass-range shown is from 0.56 to 0.58 Mq and 
the time-range is from 2.4360 x 10 6 to 2.4375 x 10 6 yr. More details of AGB evolution is given m 
the review article by Herwig jsj 
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FIG. 3: Evolution of the 12 C and 16 (in mass fraction) in the envelope (i.e., at the surface) as 
a function of total stellar mass (sequence ET13). Since the total stellar mass decreases with time 
the figure shows the time evolution of the envelope 12 C abundance. The intershell material that is 
dredged-up to the envelope contains more C than O, which eventually leads to C-star formation 
(C/O > 1). Integration of the surface abundance over the mass lost gives the 12 C yields. 
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FIG. 4: Rate selection for stellar evolution He-shell flash calculations. The numbers correspond to 
the sequence numbers in this paper. The red (dark) area corresponds to the NACRE range for the 
14 N(p, 7) rate, and the green (light) area covers the revised range for this rate, for the temperature 
range considered in this work (cf. Tabled). 
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FIG. 5: Overview of 12 C yields as a function of the 14 N(p, 7) and triple-a rate, for our M = 2 M , 
Z = 0.01 TP-AGB models as a function of nuclear reaction rates. Each point refers to the yield of 
one full stellar evolution model sequence. Lines connect points with the same triple-a rate. 
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FIG. 6: Like Fig-El for sequences with different nuclear reaction rates (Fig.0J. The yields of these 
sequence are given in TablelTTI and IIII1 
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